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The form factor provides a convenient way to describe properties of topological solitons in the 
full quantum theory, when semiclassical concepts are not applicable. It is demonstrated that the 
form factor can be calculated numerically using lattice Monte Carlo simulations. The approach is 
very general and can be applied to essentially any type of soliton. The technique is illustrated by 
calculating the kink form factor near the critical point in 1+1-dimensional scalar field theory. As 
expected from universality arguments, the result agrees with the exactly calculable scaling form 
factor of the two-dimensional Ising model. 



PACS numbers: ll.15.Ha, 11.27.+d 
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Topological solitons play an important role in a wide 
range of physical systems [iHSl , and they have been stud- 
ied extensively both experimentally and theoretically. 
Although in many of these systems quantum mechan- 
ical effects are significant, theoretical studies have been 
mainly limited to the classical limit in all but the simplest 
cases. In principle, one can use perturbation theory to 
calculate 'semiclassical' quantum corrections to classical 
quantities However, this only works when the quan- 
tum effects are small and in practice can only be used 
for very simple models Q. A fully quantum mechanical 
analysis requires not only a new calculational method but 
also a different set of observables; the simplest classical 
observables do not have well-defined quantum mechani- 
cal counterparts. 

In this paper we investigate one such observable, the 
soliton form factor. The form factor is a fully non- 
perturbative quantum observable which does not rely on 
any semiclassical concepts, and it can be defined for any 
soliton in an analogous way. Form factors are used in 
many areas of physics to characterize properties of quan- 
tum objects, from atomic Q and nuclear physics Q to 
integrable systems Q. The soliton form factor describes 
the scattering of a particle with a soliton; it can be loosely 
interpreted as (the Fourier transform of) the soliton pro- 
file in the quantum theory. It is therefore the most natu- 
ral quantum observable beyond the soliton mass, and it 
carries a large amount of non-trivial information about 
the soliton and its interactions. By studying the form fac- 
tor, one can therefore move away completely from semi- 
classical ideas of soliton shape to work with fully non- 
perturbative results for excitations and interactions. 

Choosing a concrete example with both nontrivial crit- 
ical behaviour and well-understood semiclassical limits, 
we shall focus on the kink form factor [9|-|l3| in 1-1-1- 
dimensional field theory. Previously [1J|, we studied 
semiclassical aspects of kinks with lattice Monte Carlo 
simulations by measuring the field correlation function 
in the presence of a kink. Here we show how that the 
same approach can be further developed to calculate the 
kink form factor in a fully non-perturbative way. We 
demonstrate this by carrying out simulations near the 
critical point, where we find excellent agreement with 



exact results from the two-dimensional Ising model, as 
predicted by universality arguments. This approach can 
be generalized to solitons in other theories. 

Let us consider a theory with a real scalar field 4>{t, x) 
and kinks in 1-|-1 dimensions. Our discussion will be valid 
for any such theory, but as a concrete example, we use 
the X(j)'^ model with Lagrangian 



(1) 



i(a,0)(5'^0)-f ^0^-^04. 

For fixed coupling A there is a critical mass parameter 
value below which the Z2 symmetry of this theory is 
spontaneously broken, and the scalar field has a vacuum 
expectation value {(f>) = ±w. In the classical theory = 
and V = my^G/X. The kink corresponds to a state 
which interpolates between one vacuum on one side and 
the other vacuum on the other side. 

We assume that the system is in a state with one 
kink. Classically, this simply corresponds to the exact 
kink solution 0kink(a^) = v tanh(ma;/\/2). In the quantum 
theory, the same can be achieved by imposing twisted 
antiperiodic boundary conditions in the spatial direc- 
tion [1^. In fact, this only restricts the number of kinks 
to odd values, but states with more than one kink are ex- 
ponentially suppressed in the infinite volume limit. Note 
that this way of preparing the system preserves transla- 
tion invariance and is fully non-perturbative as it makes 
no reference to any classical background configuration. 

The ground state |0) of this one-kink sector of the the- 
ory corresponds to a kink in a momentum eigenstate with 
zero momentum, and its energy Eq is therefore simply the 
kink mass M. This state is therefore translation invari- 
ant. Above this, the spectrum consists of moving kink 
states \k) with momentum k and energy 
up to the energy Eexc of the first excited state of the kink. 
Above M + m, there are also states with one or more free 
scalar particles. 

We want to calculate the kink form factor f{k,k'), 
which is defined as the matrix element 

f{k,k')^-{k'\m\k), (2) 

V 

where we have scaled it by the vacuum expectation value 
V to make it independent of the field normalisation, and 
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the momentum states | k) have the Lorentz invariant nor- 
malisation 



(fc'Ifc) =2n6{k~ k')Ek. 



(3) 



Lorentz invariance of the theory implies that, when ex- 
pressed in terms of rapidities Pk = arcsinh k/M , the form 
factor is a function of the rapidity difference only [lol |. 

f{k,k') = f{f3k-M- 

Semiclassically, the kink form factor is given by the 
Fourier transform of the static kink solution [ll| , 



fc\{f3) = -inv^ 



sinh ^TTW^/J 



(4) 



This means that even in the quantum theory the form 
factor can be thought of as the effective kink profile. 
However, this interpretation should not be taken literally 
because, as always, there are many quantum observables 
that have the same semiclassical limit. The semiclassical 
approximation is valid at weak coupling. In our model 
(P), the dimensionless coupling is X/m^ = Q/v"^, thus 
weak coupling implies a large vacuum expectation value. 

What makes the theory ^ a particularly useful test 
bed is that the form factor is also known exactly at 
strong coupling, by which we mean near the critical point 
w in the quantum theory. The theory is in 
the same universality class as the two-dimensional Ising 
model, and near the critical point the form factor should 
approach the exact Ising model result [H, [l3| , 



/ising(/3) = icoth^. 



(5) 



Matrix elements like ([2]) cannot be computed directly 
using Monte Carlo simulations. Instead, the basic ob- 
servable is the field correlation function, which we con- 
sider in the ground state |0) of the one- kink sector. We 
calculate it in momentum space, taking the Fourier trans- 
form in space but not in time, and write a spectral ex- 
pansion in terms of energy eigenstates |a) with energies 

Ea, 

(6) 

Lattice Monte Carlo simulations are carried out in Eu- 
clidean space, which is obtained by carrying out a Wick 
rotation t — > it. This does not affect the coefficients of 
the spectral expansion, but the exponentials become real. 



WO,k)<j,{t,q)) 



(O|0(fc)|c^)(a|0(g)|O) 
(0|0) 



-t(E„-Eo) 



At long enough time separation, 

1 

t > 



(7) 
(8) 



the dominant contribution comes from the single-particle 
moving kink states \k). For them, the coefficient of the 
expansion is essentially the form factor, because 



{k'\<j>{q)\k) = vf{k, k')2Ti5{k -q- k'). 



(9) 



The momentum conservation delta function restricts the 
expansion to only states with overall momentum fc, and 
therefore we have 

L EkEo 



+0 (e-*(^- 



-A-f) 



(10) 



where L is the spatial length of the system, and we have 
used (0|0) = LEq as implied by our normalisation. 

Furthermore, the Euclidean spacetime is necessarily fi- 
nite in actual Monte Carlo simulations. We assume pe- 
riodic boundary conditions in the time direction, and 
denote the length of the system by T. In the periodic 
Euclidean time, the field correlator is 

where U{t) = exp{—Ht) is the Euclidean time evolution 
operator. As in Eq. (O, at long enough time separations 
the only contribution comes from single-particle kink mo- 
mentum eigenstates \k), so we can approximate the trace 
in Eq. (1111) by an integral over them. 



(<^(0,fc)</>(i,<z)) = 



I^{k'\U{T-mq)U{mk)\k') 



I 2^ik'\U{T)\k') 



Using 



{k'\U{t)\k) = 2Tr5{k~ k')Eke 



we can write the denominator as 



(12) 
(13) 

(14) 



Inserting complete sets of momentum eigenstates, the nu- 
merator becomes 



dk' 



ik'\U{T~t)^{q)U{t)m\k') 



o ,^ f dk' V^\f(k' ^ k,k')\^ 

^2TT6{q+k) / ■ — ii-e--^'='(^-*)--^'='-fe*. 

J 2-K Ek'-kEk' 

(15) 

As illustrated in Fig. [l] this integral has a simple geo- 
metrical interpretation: The kink travels from time to 
time t at momentum k' — k, where it interacts with a (f> 
particle. This changes its momentum to k' , with which 
it moves forward in time through the periodic boundary 
back to time 0. To calculate the integrals ((Ti)) and (|15l) . 
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FIG. 1. Illustration corresponding to Eq. (jlSp . The kink is 
constrained by periodic boundary conditions in the Euclidean 
time direction, so the worldline must match up at either end 
of the lattice. The defect is not point-like, and has a form 
factor which is represented here by the interaction between 
the scalar and the defect having a finite size. 



we use the saddle point approximation. The saddle point 
ko for Eq. (|15p is found by minimising the action 



S{k') = Ek'{T -t)+ Ek'-kt - MT 



(16) 



for given t. By approximating the integral by a Gaussian 
around the saddle point, we obtain 



(0(O,fc)0(t,g)) 
_ 2n5{k + q) It i;2|/(A;o,fco-fc)| = 



M ^Ek,-kEk,S"{ko) 



(17) 



This approximation is only valid when the action is suf- 
ficiently peaked and well approximated by a Gaussian. 
The latter assertion requires 



5"(fco)2»^(*)(fco 



(18) 



This implicitly imposes a lower limit for the system size 
T in the time direction, so for higher k we need to use 
larger lattices. As usual, the lattice size also has to be 
larger than any inverse mass, including the kink mass M . 

Finally, we note that because cf) is real and the kink has 
odd parity, the form factor is odd and purely imaginary. 
Therefore we can use Eq. ([T7|) to determine it from the 
field correlator, up to a sign. For given k and i, we obtain 
the saddle point fco by minimising Eq. (jl6p . and the form 
factor for rapidity difference j3 is given by 



/(/3)=/(fco,fco-fc) 



^{4>{{),k)4>{t,-k)) ( MEk,-kEk,S"{ko)V'^ 



where 



T 



p = arcsmh— -: — arcsmh— — — — . 

2M 2M 



While Eq. p9)) is an approximation, it becomes exact 
for sufficiently large T as discussed above; one must also 
satisfy Eq. dH) by excluding small t. 

We tested this result by calculating the form factor 
near the critical point using lattice Monte Carlo simula- 
tions. The Euclidean lattice action for the theory ([1]) is 
given in lattice units by 



(19) 



(20) 



+ 2- — 



4! 



(21) 



with A = 0.6. Square lattice sizes oi L ^ T E 
{125, 250, 375} were used. A kink is created by imposing 
antiperiodic boundary conditions in the space direction, 
4){x + L,t) = —(f){x,t). This also leads to discretisation 
of momentum, k = (2n + 1)tt/L. 

We measured the momentum space unequal-time field 
correlator at various time separations and then used 
Eq. ((T9)) to calculate the form factor for various rapidities. 
At the strong couplings used here, the correlator mea- 
surements are reliable even at very long distance thanks 
to the hybrid Monte Carlo algorithm which was helpful 
in fighting critical slowing down and thermalizing long- 
distance modes efficiently. In principle, Eq. (jl9p gives the 
form factor for a range of /3 from a single choice of param- 
eters {k, L, m^} because the same simulation gives the 
correlator for all values of the time separation t. However, 
these values are strongly correlated; we report only one 
data point per combination with a quoted error obtained 
from a bootstrap resampling of all measurements [l§] . 

Given Eq. ([5]), the time separation t has to be long 
enough that excited states and two-particle states are 
suppressed sufficiently. This happens when t > 1/2M. 
However, at greater distances statistical noise starts to 
grow. Therefore, we select the value of t with the smallest 
statistical error within the permitted range. 

In addition to the field correlator, Eq. (|19l) also in- 
volves the vacuum expectation value v of (j), and the kink 
mass M. We measured v using simulations with peri- 
odic boundary conditions. We take the ensemble average 
V — {\J2xev ^^(^)|) of the average field's absolute value. 

To obtain the kink mass M, we again used Eq. and 
the saddle point approximation, but this time taking the 
lowest available momentum, k — n/L, leaving the higher 
momentum measurements as independent datasets for 
study of the form factor. Then, as long as fc ^ M, 
Eq. ((TE)) simplifies and we find the saddle point ko = kt/T 
for arbitrary t. We can, therefore, apply the saddle point 
approximation analytically, and we find [l3| 



((/.(0,fc)(/.(t,g)) oc ^-V^^^,t-VM-+ik-ko)HT-t)+MT^ 

(22) 

We then fit the k — n/L correlator data to this expression 
to obtain the mass M with a bootstrap error. The results 
of this fitting are shown in Fig. [51 Alternatively, the mass 
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FIG. 2. The kink mass M as a function of m? in different 
volumes. The data agree with the known results for the Ising 
universality class extrapolating to the infinite- volume 

limit, the dependence on is linear. 



FIG. 3. The form factor as a function of the rapidity dif- 
ference. Measurements are shown for several lattice sizes at 
various momenta, and the form factors for the Ising model 
and for a semiclassical kink are also shown. 



could also be calculated from the free energy difference 
between the kink and vacuum sectors [isj . 

This way, we have measured all the quantities that 
appear in Eq. ([TO]). We can calculate the form factor 
/(/?) at a wide range of rapidities for different momenta 
and different values of m^. Finally, we have checked the 
consistency of the saddle point approximation leading to 
Eq. dni) when Eq. (HH) is satisfied. 

The results in the critical regime are shown in Fig. |31 
together with the exact Ising model result ([5]) for com- 
parison. The agreement is very good. This demonstrates 
that we can calculate the kink form factor reliably even at 
strong coupUng where perturbative approaches fail. Un- 
Hke the Ising model, the scalar field theory ([Ij in which 
we carried out the calculation is not exactly solvable, 
and we made no use of any special features of the the- 
ory. Therefore we expect that the same method will work 
equally well in other theories. 

We have shown how to calculate soliton form factors 
nonperturbatively from field correlation functions mea- 
sured in lattice field theory simulations. The approach 
can be applied directly to other theories with kinks, such 
as the Sine-Gordon model, and generalisation to other 
theories where twisted boundary conditions can create 
topological solitons [20l - [2^ should be straightforward. 

In more complicated theories, one will obtain several 
form factors which describe interactions of the soliton 
with different particle species. In the case of point-like 
solitons, such as 't Hooft-Polyakov monopoles in 3-1-1- 
dimensional gauge field theory [g^j, the calculation will 
follow the same lines. Generalisation to extended soli- 
tons, such as domain walls, strings or higher-dimensional 
membranes is less trivial but should still be possible. 
For example, this will make it possible to study non- 



perturbatively the quantum mechanical properties and 
interactions of cosmic strings. 
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